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Abstract 

Crossing probabilities for critical 2-D percolation on large but finite lattices have 
been derived via boundary conformal field theory. These predictions agree very well 
with numerical results. However, their derivation is heuristic and there is evidence of 
additional symmetries in the problem. This contribution gives a preliminary exami- 
nation some unusual modular behavior of these quantities. In particular, the deriva- 
tives of the "horizontal" and " horizontal- vertical" crossing probabilities transform 
as a vector modular form, one component of which is an ordinary modular form and 
the other the product of a modular form with the integral of a modular form. We 
include consideration of the interplay between conformal and modular invariance 
that arises. 
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1 Introduction 



Percolation is perhaps the simplest non-trivial model in statistical mechanics. 
It is very easy to define, and exhibits a second-order phase transition between 
the percolating and non-percolating states. A broad array of techniques have 
been brought to bear on it over a period of many years. Its behavior is of 
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current interest, and has been studied via renormalization group, conformal 
field theory, Coulomb gas methods, computer simulation, as an example of 
super symmetry, and using rigorous mathematical methods. (For general re- 
views, see [1],[2], for some recent results including a list of references see [3]). 
In this contribution, we restrict ourselves to the study of crossing probabilities 
for two-dimensional systems at the percolation (phase transition) point p c . 

Although percolation is, as mentioned, arguably the simplest model that ex- 
hibits a second order phase transition, the ease of formulation of the model 
is in another sense deceptive, tending to conceal its inherent complexity. The 
wide range of approaches taken to it already attests to this subtlety. The 
ultimate reason is suspected to be the unconstrained nature of the model, 
encompassing a variety of symmetries. 

In Section 2 we review percolation and crossing probabilities. Exact analytic 
expressions for these quantities are known from boundary conformal field the- 
ory. We transform these results into a form suitable for the present analysis. 
Section 3 briefly introduces modular forms. Then the unusual modular be- 
havior of the derivatives of the crossing probabilities is examined. The results 
here are preliminary; a full treatment will appear elsewhere [4] . 



2 Crossing Probabilities 

The properties that we consider here are critical and therefore universal, the 
same for a wide variety of types of (isotropic) percolation and lattice struc- 
tures. However, for definiteness, when we are specific we will refer to bond 
percolation on a square lattice. Bond percolation is defined by placing a bond 
with (independent) probability p on each edge of the lattice. Consequently, 
there are 2^ possible bond configurations with < N B < N, where N B is 
the number of bonds in a given configuration and N is the total number of 
edges. The connected bonds in each configuration form clusters. At p c , (note 
that duality implies p c = 1/2 on a square lattice), as the lattice is taken to 
infinity, one or more infinite clusters just appears. For p < p c , there is no 
infinite cluster. 

The crossing probabilities are defined by considering a finite rectangular LxL' 
lattice as L,L' — > infinity with fixed aspect ratio r = width/height = L/U. 
(Below, we will allow the shape of the lattice to change.) Then the probability 
of a configuration connecting the left side and the right side of the rectangle is 
the horizontal crossing probability Hh(r). The probability of a configuration 
connecting all four sides is the horizontal- vertical crossing probability Uh v (r). 
These quantities are known to depend only on the aspect ratio r (for a rect- 
angle) by extensive numerical work and the hypothesis (and consequences) of 
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conformal invariance. In fact, they enjoy an even wider invariance, as discussed 
below. 

Next, to motivate the conformal approach to this problem, we consider the 
Q-state Potts model. This generalization of the Ising model employs a spin 
variable Sj with values Sj = 1, 2, Q (Q = 2 corresponds to the Ising model) 
on each site i of the lattice. The (reduced) Hamiltonian is 

H = K £ 5^ (1) 

<ij> 

where < ij > denotes nearest neighbor sites. By introducing the variable 
x = e K — 1 one may rewrite the partition function as follows: 

z =y. n a (2) 

{si} <i,j> 

Expanding the product, one may perform the sum over Sj in each term. Rep- 
resenting the presence of a factor i by a bond then gives rise to a graphical 
representation of Z, known as the random cluster or bond-correlated Potts 
represent at ion [5] , [6- 9] 

Z = J2 Q Nc x Nb (3) 

G 

where the sum is over all possible graphs consisting of N b bonds arranged in 
N c clusters (counting single isolated sites as clusters). This model is known 
to have a critical point for all < Q < 4. (On the square lattice, by duality, 
%c = Q 1 ^ 2 )- For Q — 1, the set of configurations included is exactly that, 
including the weighting, as for bond percolation, with p = xj[x + 1). For 
other Q values, the configurations are the same but weighted differently. In 
addition, Eq. (3) allows us to extend the number of states to Q G R. Thus we 
can envision a continuous change from the Ising model, say, to percolation. 
Further, the central charge is known as a function of Q; in particular, c = 
for Q = 1 (critical percolation). 

In order to study the crossing probabilities, following [10], consider for def- 
initeness 11^ on a rectangle. Let be the partition function of the Potts 
model with the spins on the left vertical side fixed in state a and those on the 
right vertical side fixed in state b. The spins on the rest of the boundary are 
unrestricted. Then 

q-A zj 
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(4) 



with a ^ b. Note that this expression differs from the one in [10] due to the 
normalization of Z. The expression makes no sense for Q = 1, of course, but 
it allows a solution to the problem if we first express the partition functions 
using boundary operators and then take the limit. 

Since we are at criticality, the partition function can (in the limit of a large 
lattice) be expressed using conformal field theory. The critical Potts models 
are known to correspond to a certain series of minimal models. A change of 
boundary conditions is introduced by means of a boundary operator [11]. If 
the coordinates of the corners of the rectangle are zi, z 2 , z 3 , z 4 we thus have 

Z ab = Zf {4>fa{zi)(t>af{z2)(t>fb{z-i)4>bf{Zi)) (5) 



where Zf is the partition function with free boundary conditions and the 0s 
are boundary operators. The next step is to identify <p a f. This is done by 
comparison with known results for the Ising and Q = 3 state Potts models. In 
these cases, the operator that changes between fixed boundary conditions a 
and b is known to be 0ri,3)- On the other hand, one can implement this change 
by bringing together two points Z\ and z 2 where the boundary conditions go 
from a to / and / to b, respectively. Using the operator product expansion then 
gives a term that must be the operator in question. By the fusion rules, the 
only operator that can satisfy this is seen to be <p a f = (f>n t2 ) ■ This argument 
is doubly satisfying, since the conformal dimension of </>(i,2) is h — in the 
limit Q — > 1, which is a necessary requirement for the crossing probability 
to be conformally invariant. Further, the operator is level two, so that the 
differential equation satisfied by its four-point function is second order. 

It is conventional to consider the problem on the upper half plane, which 
may subsequently be mapped onto a rectangle via the Schwarz-Christoffel 
transformation, with the four corner points taken as images of — 1/k, — 1, 1,1/k 
(with < k < 1). Then the crossing is between the intervals —1/k < x < —1 
and 1 < x < 1/k on the real axis. The four-point correlation functions depend 
only on the cross-ratio 

= (X4 - x 3 )(x 2 - x{) = ( 1 - £: \ 2 . . 

(x 3 - x l )(x 4 - x 2 ) \l + k) 



One then finds, by means of standard conformal manipulations, that the cor- 
relation function satisfies a Riemann equation with the two solutions F(X) = 1 
and F(X) = A 1/3 2 i 71 i(l/3, 2/3; 4/3; A). One can pick the correct linear combi- 
nation by imposing the physical constraints that F-^OasA^O(r^oo) 
and F — >■ 1 as A — > 1 (r — ► 0). The result is 

Uh{x) = mw xl/32Fi{1/3 ' 2/3; 4/3; A) (7) 
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Mapping Eq. (7) onto a rectangle, one finds that the aspect ratio becomes 
r = 2K/K', where K' and K are the complete elliptic integrals. This result 
has been extensively tested via Monte Carlo simulations [12], [13], and there is 
little doubt as to its correctness. Note that one can transform it to any compact 
shape with four identified points, not just a rectangle; the crossing probability 
will be the same as the crossing on the half-plane with the corresponding 
half-plane cross-ratio. Thus, for instance, by consideration of the Schwarz- 
Christoffel formula it is easy to see that the crossing probability on a rhombus 
of any angle is the same as on a square. This point has been investigated 
numerically in Fig. 4.1 of [13]. Of course the same invariance holds for the 
"horizontal- vertical" crossing considered below, only the function F changes. 

In general, when one makes a conformal transformation z — > w(z) of a cor- 
relation function, factors of (w'(z)) h appear. In addition, transforming from 
the upper half plane to a shape with corners, the partition function gains a 
(non-scale invariant) factor L ac , where L is the length scale, c is the central 
charge and a depends on the geometry. Similarly, a correlation function with 
boundary operators sitting at the corner points gains a factor L~^^ h , where 
7 is the interior angle at the corner [14]. The last two (non-scale invariant) 
effects occur because the transformation is only piecewise analytic and has 
singular points at the corners. However, for critical percolation, c = = h, 
and all of these factors become unity. Thus the crossing probabilities are in- 
variant under transformations that are conformally invariant in the interior of 
a region and only piecewise conformally invariant on its edges. 

The "horizontal- vertical" crossing probability Hhv m ay be obtained similarly 
[15]. We omit most details of the argument. A four-point boundary operator 
correlation function arises once again. The main differences with 11^ are the 
complexity of the boundary conditions and the absence of any direct identifi- 
cation of the boundary operator. Instead, one considers low-lying null vectors 
in the c = = h Verma module. Through level 5, there is only one which 
leads to solutions that satisfy the physical requirements. These are 



n ftw (r) 
n ftw (r) 



n ftw (i/r) 
n fc (r) 



(8) 



which translate into 



F(X) 
F(X) 



F(l-A) 

n h (A) 



(9) 
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respectively. Applying these conditions to the level 5 solutions, one finds 

IMA) = AV3 2^i(l/3, 2/3; 4/3; A) - t^A 3 F 2 (1, 1,4/3; 2, 5/3; A)(10) 



where 3F2 is a generalized hypergeometric function. The first term is just 
Uh and the second subtracts configurations with horizontal but no vertical 
crossings. 

The differential equation satisfied by F may be written 



d\ 2 



(A(A-1)) + 



1 d 
2\-ld\ 



(2A - 1) 



F = 
(11) 



The factorized form exhibited in the second line is of interest since 1, n^, and 
Uhv span the solutions of the equation formed by letting the rightmost factor 
act on F alone, i.e. 



(12) 



In what follows, it is convenient to consider the r-derivatives of Il h and H hv 
(on the rectangle) which we will denote U.' h (r) and H' hv (r). Note that H' h , 
for instance, is interpretable physically as the probability density that the 
maximum horizontal extent of a cluster attached to one vertical side of an 
infinitely wide rectangle of unit height is greater than r [16]. Additionally, 
since the r-derivative is proportional to the A-derivative, Eq. (12) reduces to 
second order. 

We next proceed to express U' h and U' hv on the rectangle as functions of r, 
using the result for the cross-ratio [16, 17] 



(13) 



where $2 and $3 are the elliptic theta-functions and q = e wr (note that 
q is the square root of the usual q). Eq. (13) follows by applying Landen's 
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transformation to r = 2K/K' and Eq. (6), resulting in 



r= *(^ (14) 



where K is the complete elliptic integral written as a function of the modulus. 
Expressing the latter in terms of a ratio of theta-functions ([18] 8.197.1,2) one 
obtains Eq. (13). 

We also note, for future reference, some identities involving A. One has A = 
le ^gy 1 ' , 1 - A = and A' = 16^0^, where rj is the 

Dedekind 77-function (see Eq. (18)) with q = e 2ntT , and the differentiation is 
with respect to the independent variable r = ir. 

Eq. (13) makes it possible to re-write the differential Eq. (12) for F directly 
in terms of the aspect ratio r. One obtains 

g + a(r)! + &(r)/ = (15) 
where /(r) = F'(r), the r-derivative of the four-point function, and 



3A" 5A' 5A' 
° (r) = + 3(1^1) + 3A 

5A" /A"\ 2 A w 5A" 4(A') 2 
6(r) = -3A +3 Uj -^-^T) + 3T(X3T) ( 16 ) 



where differentiation is with respect to r. One may identify two independent 
solutions of Eq. (15) as 



^ 4 



h= v(q 2 ) 



h = \\m)f - fw (17) 



where 77 is the Dedekind 77-function 

00 

v(q) = HO- - 9 n ) ( 18 ) 

71=1 
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and fw is an even, and apparently new, function of q. Its first few terms are 
given by 




13568 
4301 




458070 
124729 




(19) 



The solution fi is seen to be proportional to Ii' h by the above, or directly [17]. 
We note for completeness that 



The connection with Tl' hv involves including a term proportional to f 2 . This 
is specified below. Note that, because of the way we have defined the aspect 
ratio r, our Il h coincides with II„ in [16, 17]. 



3 Modular Forms 

This section follows the excellent introduction given in [19]. 

A modular function or form assigns a complex number G(A) to each lattice 
A. Here the term 'lattice', following mathematical usage, refers to an infinite 
regular array of points, defined by the basis {001,002} A = Z001 + Zoo 2 . In 
addition, for any complex number A 7^ 0, G satisfies G(A) = X k G(XA), where k 
is some integer, called the weight (with k = for modular functions). Because 
of this, on dividing by oo 2 , we see that G is completely specified by g(r) = 
G(Zr + Z), where r = 001/002. Since g is also even in r, we may restrict r 
to the upper half plane. In addition, there are analyticity and growth (as 
r — > ioo) conditions on g. 

The modular properties arise on considering a change of basis. We can replace 
{001, oo 2 } by {lo[, lo' 2 } = {auoi + boo 2 , cuo\ + doo 2 } with a, b,c,d G Z, ad — bc = ±1 
without changing the lattice. This implies that g must satisfy the modular 
transformation property 



Restricting r to the upper half plane imposes ad — be = +1. The group of 
matrices implementing such transformations is the (full) modular group r ± . 




2 7 / 3 7T 2 



(20) 



>/3r(i/3) 3 




(21) 
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The matrix T 
( 



'11 



V 



-1 

1 



1 

for t — > — 1/r, generate IY 



, which implements r — > r + 1, together with S 



It follows from Eq. (21) that modular forms of a given weight are a vector 
space over C. The dimension of this space is very small, in general, a fact 
which leads to some very non-trivial relations. 

In what follows, we examine the modular properties of fi and f 2 , considered as 
functions of r = ir, with r complex. (Here it is useful to envision the rectangle 
rotated by 90°, since r is in the upper half plane for real r.) In order to 
make sense of this, we need to understand the physical meaning of a modular 
transformation in the context of percolation crossing probabilities. We do not 
include a full treatment of this question here, but restrict ourselves to the 
case of a conformally invariant physical quantity II, which is also a modular 
function (i.e. invariant under any 7 G 1^). Suppose further that II is initially 
defined on a rectangle of aspect ratio r. Now the modular transformation, as 
indicated, acts only on the basis vectors {u>i,ui 2 }. (For a rectangle, RejcJi} = 
= Im {^2})- To construct a conformal map of the rectangle, we consider 
these as displacements from a fixed origin at 0. Thus the map must satisfy 
{0,cji,cj 2 } {O,^,^} = {0,a&i+bui 2 ,cui 1 +du! 2 }. This can be implemented, 
for instance, by a projective transformation 



az 

w 



ez + 5' 
a = (uj 2 - WiViu4 



e = uj' x uj 2 — WiOJ 2 , 



5 = ijjxujiiu'i — oj'i) (22) 



which will take the rectangle into a figure with sides that are straight or arcs 
of circles. Note that the necessary conditions uo\ 7^ uj 2 and uj[ 7^ uj' 2 imply 
a, 5 7^ 0. This choice of map also has the advantage of always being 1-to-l. 
In addition, one can show that it preserves the group structure. However, the 
map is explicitly dependent on the basis vectors. Note that the quantity n 
will remain invariant. The parameter r, defined as a ratio of displacements, 
transforms in the usual way, so that n(7(r)) = n(r). Thus for each 7, modular 
invariance corresponds to conformal invariance under a particular map. Of 
course we can also begin with the physical quantity defined on a parallelogram, 
removing the conditions on {oji,cj 2 }- Note, however, that the derivative of 
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w at either of the three corners differs from the derivative of the modular 
transformation. 



We now return to the crossing probabilities. We focus on the functions fi 
and f 2 , solutions of Eq. (15), considering them as functions of r = ir, so 
that q = e mT . We have already specified the connection between fi and YY h . 
As mentioned, the solutions of Eq. (12) (and therefore those of Eq. (15)) 
span U' h and H' hv . Thus it only remains to find the correct linear combination 
of fi and f 2 . Now IY^, must, by the first of Eqs. (9), be invariant under 
W : /(t) -> T- 2 f(-l/r). The linear combination / 2 - f /i with C = 
satisfies this condition. Combining this with Eq. (20) we find that 



2 7/3 



7T 



2 



24 



n - w = -Tm^r [ " w )J + vs m (23) 



Now this function behaves like a weight 2 modular form under the operations 
t — > — 1/r and r — > r + 6, but these transformations generate a subgroup of 
infinite index of SX 2 (Z), and 11^ is not a modular form. However, the vector 

transforms under T 2 and S, which generate the theta-group T g , 



F = 



fi 



according to 



F(r + 2) 



uj 
1 



F(r) 



r- 2 F(-l/r) 



-1 
-C 1 



F(r) 



(24) 



where C is as above and u = e 2m / 3 . Thus fx is a weight 2 modular form (with 
multiplier system, cf. [19]) on T e , but f 2 is a kind of "second order" modular 
form, i.e. instead of 71/2-/2 = for any 7 e r#, (where | denotes the modular 
weight 2 operation of 2 x 2 matrices) we have that 7I/2 — f 2 is a multiple of 
/i- 

Another way to look at this situation is to consider the differential equation 
(15). For convenience we take the independent variable to be r = ir. Then, 
writing / 2 = uf\, we find a second-order equation fiu" + (2f[ + aj\)u' = for 
u'. Now a = -h'/h, where h = A' 3 A~ 5/3 (1- A)~ 5/3 so that «"/«' = -a-2f'J ' f\. 
Hence u' = 2/3h /f 2 . Using Eq. (17) for /1 and the identities for A above, we 
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have finally 



2 t)(t/2) 8 t](2t) 

3 r](T) 12 
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(25) 



It follows from the modular properties of r\ that u' is a modular form of 
weight 2 on Yg, so that /2 is the product of a modular form and the integral 
of a modular form. 

Our examination of the modular properties of j\ and ji reveals a close connec- 
tion between the two functions. Now modular forms, as mentioned, are vector 
spaces of very small dimension. Thus it is likely that one can derive f'2 given 
/1 and some physical conditions. This would be very interesting. Whether it 
is indeed possible remains to be seen, however. We have succeeded in showing 
that f\ follows (up to an overall multiplicative constant) if one assumes that 
only one conformal block contributes to it [4] . 

One can also show that the function u corresponds to the Weierstrass (- 
function on the elliptic curve Y 2 = X 3 — 1728. 
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